Draft version May 16, 2008 

Preprint typeset using I^T^K style cmulateapj v. 10/09/06 



oo 
O 
O 
(N 



6 



ENERGETICS OF X-RAY CAVITIES AND RADIO LOBES IN GALAXY CLUSTERS 

William G. Mathews 1 and Fabrizio Brighenti 1,2 
Draft version May 16, 2008 

ABSTRACT 

We describe the formation and evolution of X-ray cavities in the hot gas of galaxy clusters. The 
cavities are formed only with relativistic cosmic rays that eventually diffuse into the surrounding gas. 
We explore the evolution of cavities formed with a wide range of diffusion rates but which are otherwise 
similar. In previous numerical simulations, in which cavities are formed by injecting ultra-hot but 
non-relativistic gas, cavity formation contributes thermal energy that may offset radiative losses in the 
gas, thereby helping to solve the cooling flow problem. Contrary to these results, we find that X-ray 
cavities formed solely from cosmic rays have a global cooling effect. Most cosmic rays in our cavity 
evolutions do not move beyond the cooling radius and, as the cluster gas is displaced, contribute to 
a global expansion of the cluster gas. As cosmic rays diffuse away from the cavities, the nearby gas 
becomes buoyant, resulting in a significant outward mass transfer within the cooling radius, carrying 
relatively low entropy gas containing cosmic rays to outer regions in the cluster where it remains for 
times exceeding the local cooling time in the hot gas. This post-cavity mass outflow due to cosmic 
ray buoyancy may contribute significantly toward solving the cooling flow problem. For example the 
mass inflow in the Virgo cluster due to radiative cooling can be balanced by buoyant outflow if only 
a fraction ~ 0.0005 of the accretion energy onto the central black hole inflates X-ray cavities with 
cosmic rays. We describe the energetics, size, stability and buoyant rise of X-ray cavities in detail, 
showing how each depends on the rate of cosmic ray diffusion. 
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1. INTRODUCTION 

The hot gas in galaxy clusters loses energy by X-ray 
emission but does not cool to low temperatures. In re- 
cent attempts to solve this cooling flow problem it has 
been almost exclusively assumed that feedback energy 
from accretion onto cluster-centered massive black holes 
can be delivered to the hot gas in a manner that main- 
tains the observed temperature and density profiles in 
spite of radiation losses. The solution to this problem 
is not straightforward. Even when heat is supplied to 
the cluster gas in an ad hoc, optimized, ideally fine- 
tuned manner, either concentrated or distributed over 
many tens of kiloparsecs, it is found that the temperature 
and density profiles deviate strongly from those observed 
(Brighenti & Mathews 2002; 2003). For example, highly 
idealized flows in which radiative cooling is perfectly bal- 
anced by local heating at every radius are inconsistent 
with the secular increase in gas density associated with 
stellar mass loss in the cluster-centered galaxy. In most 
galaxy groups and clusters a positive temperature gradi- 
ent (dT/dr > 0) is observed in the inner regions. Since 
the coolest gas is closest to the central source of AGN 
heating, this central gas must be heated with exquisite 
precision (on short timescales) to maintain the low tem- 
peratures observed. 

In spite of these difficulties, a variety of heating mech- 
anisms continue to be investigated to understand the en- 
ergetics of hot cluster gas. Heating is usually assumed to 
be associated with the formation of X-ray cavities by jet 
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feedback energy proceeding from the central black hole. 
Weak shock waves emerge when the cluster gas is dis- 
placed as cavities form (see McNamara & Nulsen 2007 for 
a review). Weak shocks, observed in a few clusters (e.g. 
Perseus: Fabian et al 2003; Virgo: Forman et al. 2005; 
2007), have the desirable capability of dissipating AGN 
feedback energy over large regions of the cluster g as , as 
emphasized by Fabian et al (2003). However, energy dis- 
sipation by outwardly propagating shock or sound waves 
is disproportionally concentrated in the central regions 
of clusters where the gas density gradient is smallest. 
Over time the cumulative wave dissipation in this cen- 
tral gas (through which all the waves must pass) causes 
its temperature to become hotter than observed in any 
cluster (Fujita & Suzuki 2005; Mathews, Faltcnbacher & 
Brighenti 2006). It may be possible to forestall cooling 
in a limited region of a particular cluster with a proper 
choice of shock Mach numbers (or sound wave frequen- 
cies and amplitudes), but this level of fine-tuning seems 
contrived. 

X-ray cavities are thought to provide a convenient mea- 
sure of the amount of feedback power delivered to the hot 
gas. Cavities are observed in ~ 20 — 25% of X-ray bright 
clusters (Birzan et al. 2004; Raffcrty et al. 2006), but the 
incidence of cavities increases as the cooling time of the 
central gas decreases (Dunn & Fabian 2006). The work 
done in displacing a volume V of cluster gas at pressure 
P is PV and the energy of the material inside the cavity 
is PV/(7 — 1). This corresponds to a total cavity energy 
of E cav = [7/(7 - \)]PV = APV where 7 = 4/3 is often 
chosen because many (young) cavities are filled with ra- 
dio synchrotron emission from relativistic electrons and 
there is little or no observable evidence for hot thermal 
gas inside cavities (e.g. McNamara & Nulsen 2007). 
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Assuming that the local buoyancy time t ouoy in the 
cluster gas is a measure of the cavity lifetime, Raffcrty 
et al. (2006) estimate the "cavity jet power" 

E 

*cav — T • 
^buoy 

Raffcrty et al. find that P ca v exceeds the total X-ray lu- 
minosity L x within the cooling radius (defined where the 
cooling time is 7.7 Gyrs) in about 60% of their sample 
clusters with cavities. They suggest that clusters with 
P C av < L x - and the many clusters with no known cavi- 
ties - are in a phase of their feedback energy cycle that is 
just now recovering from a recent heating episode. Raf- 
ferty et al. suggest that all the energy within the cavities 
PV/(l ~ 1) ma y be available to heat the cluster gas, not 
just the PV work done by the expanding cavity. This 
implies that the (thermal) energy of material within the 
cavity is ultimately shared with the cluster gas. Fur- 
thermore, for those clusters in which the total cavity 
jet power is favorable (P cav > L x ), it is assumed that 
this energy is distributed throughout the cluster gas in a 
manner that preserves the characteristic gas temperature 
and density profiles observed in galaxy clusters. Implicit 
in this assumed energy distribution is the requirement 
that the cavity-forming energy is delivered to the cluster 
gas in an approximately isotropic manner relative to the 
central black hole. 

Another less often considered possible resolution to 
the cooling flow problem is the hypothesis that gas is 
only heated at or near the central black hole and then 
is buoyantly transported far out into the cluster gas, i.e. 
a circulation flow (Mathews et al. 2003; 2004). One 
of the initial motivations for mass circulation was the 
two-temperature (and therefore buoyant) flow observed 
in the galaxy group NGC 5044 (Buote et al. 2003). How- 
ever, Temi, Brighenti & Mathews (2007a, b) have recently 
found far-infrared emission from dust extending out to 
5-10 kpc around many group-centered, X-ray luminous 
elliptical galaxies. Since the dust lifetime to sputtering 
destruction in the hot gas is only ~ 10 7 yrs, this obser- 
vation provides additional strong evidence for ongoing 
buoyant outflow from dust-rich cores in the central galax- 
ies. Additional support for buoyant outflow is provided 
by the large regions of iron enrichment in the hot gas sur- 
rounding cluster-centered elliptical galaxies (De Grandi 
et al. 2004). Since the 50-100 kpc size of these regions 
enriched by Type la supernovae greatly exceeds that of 
the central galaxy where they occur, an outward mass 
transfer is essential for their formation. Buoyant outflow 
is desirable since it can also preserve the observed gas 
density and temperature profiles (Mathews et al. 2004). 
Alternatively, these profiles can be preserved if mass cir- 
culates outward in momentum-driven, mass-carrying jets 
(Brighenti & Mathews 2006). 

In most theoretical studies of X-ray cavity evolution, 
it is assumed that the cavities are inflated with ultra- 
hot gas (occasionally but not always with 7 = 4/3) (e.g. 
Briiggen & Kaiser 2002; Reynolds et al. 2005; Gar- 
dini 2007; Pavlovski et al. 2007). If this heated gas 
is transported in jets from the central black hole, then 
this type of solution implicitly requires an outward mass 
circulation which in some simulations can be quite large. 
For example, in the recent double jet 3D calculation of 
Heinz et al. (2006), jets with power Wj = 10 46 ergs 



s 1 enter the cluster gas from the origin with velocity 
Vj = 30, 000 km s _1 . This translates into a mass flux of 
Mj = 2W/v 2 = 35 Mq yr^ 1 and a total injected mass of 
Mj = 3.5 x 10 9 Mq over their 10 8 year computation. If 
ambient cluster gas is entrained in the jets, the outflow- 
ing mass could be increased further. 

The jet mass flux may be very much less if most of their 
energy density is in the form of cosmic rays. While it is 
possible that jets do transport substantial masses of gas 
from the center (as in Brighenti & Mathews 2006), we 
explore here the energetics in the limiting case in which 
X-ray cavities are formed only with relativistic cosmic 
rays that can also diffuse into the cluster gas. For sim- 
plicity we assume that most of the cosmic rays in cavities 
arrive in jets and are not produced in local shocks. In 
this limit the cavities are formed with pure energy with 
no appreciable component of rest mass or momentum as 
in previous numerical simulations. Because of their dif- 
fusive nature, cosmic rays can eventually penetrate the 
cluster gas and contribute to the local pressure support. 
As a consequence the local cluster gas density is reduced 
by cosmic rays and becomes buoyant. Inhomogeneously 
distributed cosmic rays are a natural driver of buoyant 
mass outflow. Gas carried out by cosmic ray buoyancy 
may never return to the cluster center. 

However, we find that few cosmic rays move beyond 
the cooling radius, so the cluster gas expands and cools 
globally. Although some heating is expected from shock 
waves, this heating is offset by the global cooling. Con- 
sequently, X-ray cavities containing cosmic rays result in 
a net cooling of the cluster gas, not heating as generally 
assumed. Buoyant mass outflow resulting from inhomo- 
geneous cosmic rays and global expansion of the cluster 
gas may help to resolve the cooling flow puzzle. 

Although it is generally believed that X-ray cavities 
are formed with cosmic rays, this conjecture has only re- 
cently been tested with detailed calculations (Mathews & 
Brighenti 2007a, b). However, spherical steady state clus- 
ter flows including cosmic rays (e.g. Bohringer & Mor- 
fill 1988; Loewenstein et al. 1991) have been developed 
to explore the possible large scale dynamical influence of 
cosmic ray pressure gradients and the dissipation of their 
energy into the thermal gas. 

For simplicity in this initial treatment of the global 
timc-clepcndcnt energetics of cluster cavity formation 
with cosmic rays, we do not specify the physical nature 
of the relativistic fluid (electrons or protons) nor do we 
calculate radio or inverse Compton fluxes/luminosities 
which would require additional assumed parameters - 
these important details will be considered in future pa- 
pers. Nevertheless, the spatial distribution of cosmic 
rays we calculate defines the region where radio emission 
could be expected. Since we do not consider Coulomb 
heating, the cluster gas responds adiabatically to cosmic 
ray pressure gradients. 

2. EQUATIONS AND COMPUTATIONAL PROCEDURE 

The combined Eulerian evolution of (relativistic) cos- 
mic rays (CRs) and thermal gas can be described with 
the following four equations: 

§?+V-pu = (1) 
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where we suppress artificial viscosity terms. Pressures 
and thermal energy densities in the plasma and cosmic 
rays are related respectively by P = (7 — l)e and P c — 
(7c — l)e c where we assume 7 = 5/3 and j c = 4/3. The 
cosmic ray dynamics are described by e c , the integrated 
energy density over the cosmic ray energy or momentum 
distribution, e c oc / EN(E)dE oc / p 4 f(p){l+p 2 )- 1 / 2 dp. 

The first three equations are the usual equations for 
conservation of mass, momentum and thermal energy 
in the hot thermal cluster gas. We do not include op- 
tically thin radiative losses since our intention here is 
to study the energetics of cavity creation alone without 
the complicating effects of a secular radiative energy loss 
and central cooling. Note that the CR pressure gradient 
in equation 2 contributes to the motion of the thermal 
gas. This exchange of momentum between CRs and gas 
arises as the CRs diffuse through magnetic irregularities 
(Alfven waves) that are nearly frozen into the hot ther- 
mal gas. However, magnetic terms do not explicitly enter 
in the equations because typical magnetic fields in cluster 
gas ~ 1 — 10/iG (Govoni & Feretti 2004) are too small, 
i.e. the magnetic energy densities ~ B 2 /8n <> 10~ n erg 
cm~ 3 are generally much less than the thermal energy 
density in the hot gas. In addition, the Alfven velocity 

v A = B/iAirp) 1 / 2 = 2rie 1/2 B(pG) km s" 1 is typically 
much less than the sound or flow speeds in cluster gas 
so the Alfven velocity of the magnetic scatterers can be 
ignored (e.g. Drury & Falle 1986, Jones & Kang, 1990). 

Equation 4 above describes both the advection of CRs 
with the gas and their diffusion through the gas. A mass 
conservation equation for the CRs is unnecessary because 
of their negligible rest mass. The CR diffusion coefficient 
k is difficult or impossible to calculate in the absence 
of detailed information about the magnetic field topol- 
ogy which is currently unknown. However, we expect 
that k may vary inversely with the density of the ther- 
mal gas, assuming that the magnetic field strength also 
scales with density. For simplicity we ignore for now 
any dependence of n on CR particle momentum. Since 
observed radio lobes are very approximately spherical, 
we assume that k is isotropic, consistent with a highly 
irregular global magnetic field. For these preliminary 
calculations it is not necessary to specify the CR compo- 
sition, either electrons or protons can dominate as long 
as they are relativistic. Finally, we assume that the total 
CR energy density is not substantially reduced by losses 
due to synchrotron emission or interactions with ambient 
photons or thermal particles during the cavity evolution 
time. Dissipation of cosmic ray energy into the thermal 
gas is probably not important in X-ray cavities since the 
gas temperature in the cavity rims is in fact relatively 
cooler. 

The set of equations above are solved in (r, z) cylin- 
drical coordinates using a 2D code very similar to ZEUS 
2D (Stone & Norman 1992). To be specific, we study 
the evolution of X-ray cavities in the well-observed Virgo 



cluster using the analytic fits to the observed gas temper- 
ature and density profiles suggested by Ghizzardi et al. 
(2004). The computational grid consists of 100 equally 
spaced zones in both coordinates out to 50 kpc plus an 
additional 100 zones in both coordinates that increase 
in size logarithmically out to ~ 1 Mpc. We adopt re- 
flection boundary conditions for gas and cosmic rays at 
the origin and, in view of the large computational do- 
main, also at the outer grid boundary. We assume a 
spherical gravitational field g = (g z ,9r) that establishes 
exact initial hydrostatic equilibrium for the Virgo cluster 
gas pressure gradient. The cosmic ray diffusion term in 
equation 4 is solved using implicit Crank-Nicolson dif- 
ferencing. While this differencing scheme is stable for 
all time steps, we restrict each time step by the stabil- 
ity condition for explicitly-differenced diffusion as well as 
the Courant condition for numerical stability. Shocks are 
treated with a standard artificial viscosity. 

We assume that the X-ray cavity is formed by CRs that 
propagate in a non-thermal jet from the central black 
hole (AGN) to some fixed radius. The CRs are deposited 
in a gaussian-shaped sphere of characteristic radius r s — 
2 kpc located at r cav = (r, z) = (0, 10 kpc), i.e. 10 kpc 
along the z-axis. The CR source term in equation 4 is 
therefore 
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when t < t cav . The integral of (r s 7r 
space is unity. 

In the calculations described below our principal ob- 
jective is to explore the unique diffusive effects of cos- 
mic rays on cavity formation and energetics. We there- 
fore restrict the total CR energy of all calculations to 
E C av — 1 x 10 58 ergs and adopt t cav = 2 x 10 7 yrs as 
the CR injection time. We choose t cav to be consistent 
with X-ray cavity observations; it is shorter than the lo- 
cal buoyancy time in the cluster gas but sufficiently long 
not to produce strong shocks which are not commonly 
observed. At times t > t cav when S c = 0, the total CR 
energy E cr = J e c dV over the grid volume remains very 
nearly constant but changes slightly due to advection in 
adiabatic compressions or rarefactions. 

Since the CR diffusion coefficient is poorly known at 
present - and may vary from one AGN event to another 
- we consider a wide range of density-dependent coeffi- 
cients, 



10 30 crn^- 1 



n e < n e o cm 3 



10 30 (n e0 /n e ) cm 2 s 1 : n e > n e0 cm 

In general, for reasons discussed above, we assume that 
the CR diffusion k varies inversely with cluster gas den- 
sity. But k must be sufficiently large so that the CR 
density is approximately uniform in the X-ray cavities 
where the cluster gas density is lowest, i.e. we assume 
that CRs also diffuse inside the cavity. For example when 
E C av = 1 x 10 58 ergs we find that the maximum cavity 
radius is ~ 5 kpc so we require that n in the cavities 
be at least (5 kpc) 2 /t cav w 4 x 10 29 cm 2 s~ 1 and this 
condition is ensured by our adopted n(n e ) above. In re- 
gions of higher n e the density parameter n e0 determines 
the CR diffusion coefficient. In the following we consider 
n e0 = 6x 10~ 3 -6 x 10~ 6 cnr 3 . The largest CR diffusion 
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coefficient k (corresponding to n e o = 6 x 10~ 3 ) is sim- 
ilar to that required by Mathews & Brighenti (2007b) 
to explain a common age (10 8 yrs) for the large radio 
lobes and the cavity jet (thermal filament) observed by 
Forman et al. (2007) in the Virgo cluster. 

3. CAVITY EVOLUTION WITH DIFFERENT 
COSMIC RAY DIFFUSIVITIES 

Figure 1 shows the gas density p(r, z) and cosmic ray 
energy density e c (r, z) at six times during the evolution 
of an X-ray cavity and its cosmic ray (radio) lobe. In 
Figure 1 n e o = 6 x 10~ 3 cm~ 3 so the diffusion coeffi- 
cient is rather large, but identical to that used by Math- 
ews & Brighenti (2007b) to describe the evolution of the 
cavity jet (thermal filament) in M87/Virgo. The cavity 
in Figure 1 is formed 10 kpc along the (horizontal) z- 
axis during time t cav = 0.02 Gyr. At time 0.024 Gyr, 
shortly following t cav , the cosmic rays (dotted contours) 
are tightly confined inside the cavity. By time 0.066 Gyr 
the cosmic rays have diffused through the cavity walls 
forming a small radio lobe and the cavity has just disap- 
peared - we define the cavity as that region where the 
plasma density is lower than the original Virgo density 
by at least a factor 1/3. A small vortex at (r, z) — (5, 20) 
kpc is visible at this time. During the first three times 
shown in Figure 1 a (relatively cooler) thermal feature 
(the "cavity jet") is seen to rise along the z-axis. The 
evolution of this feature is described in detail in Mathews 
& Brighenti (2007b) and we will not repeat that discus- 
sion again here. During the three later times in Figure 1 
the denser parts of the cavity jet have fallen back toward 
the center of Virgo and the hot gas dynamics become 
more quiescent. However, the cosmic rays are seen to 
diffuse into a progressively larger region elongated along 
the z-axis. Although not apparent in the gas density con- 
tours, there is a net outward buoyant migration of the 
gas in the region occupied by the cosmic rays (the radio 
lobe). 

Figure 2 shows the cavity evolution when the cosmic 
ray diffusion is 1000 times lower (when n e > n e0 — 
6 x 10~ 6 cm~ 3 ) but with all other parameters unchanged. 
Now the cosmic rays are seen to be very tightly confined 
to the cavity region until at least time 0.1 Gyr when 
the cavity is still visible. In this evolution the outer 
parts of the cavity break away forming a vortex that 
migrates away from the z-axis, carrying its own cosmic 
rays. By time 0.3 Gyr the brightest parts of the radio 
image (corresponding to the largest cosmic ray energy 
density) should consist of two separate regions, a feature 
along the z-axis extending out to 40 kpc and the vor- 
tex at (r, z) = (14, 38) kpc. At this time (0.3 Gyr) we 
see an enhanced gas density that accompanies the cos- 
mic rays along the z-axis. At still later times the cosmic 
rays continue to reside mostly along the z-axis but by 
time 0.9 Gyr a region of low cosmic ray energy density 
becomes visible in the r-direction along the trajectory of 
the receding vortex (now at (r, z) — (25, 28) kpc). 

When X-ray cavities are formed in an atmosphere ini- 
tially at rest, as we assume here, the vortex region at time 
0.3 Gyr in Figure 2 would appear as a ring when viewed 
along the z-axis. To our knowledge no ring-shaped radio 
features have been observed. This could simply be due 
to the faintness of such regions since radio-synchrotron 
electrons tend to not to produce observable GHz emission 



after about 10 s yrs. Nevertheless, it is clear from Figures 
1 and 2 that the radio lobe morphology can in princi- 
ple provide valuable information concerning the cosmic 
ray diffusion coefficient about which very little is known 
at present. Finally, since gas phase metal abundances 
tend to increase toward the centers of clusters, regions 
of enhanced abundance are expected to accompany the 
outward buoyant migration of cosmic rays. 

As a further aid in interpreting Figures 1 and 2, Fig- 
ure 3 shows the pressure profiles along the z-axis at four 
times for both evolutionary calculations. At early times 
the cavity is visible as the region where P c > P. Within 
the cavity the pressure gradient (~ dP c /dr) is nearly flat 
because of the enormous pressure scale height of the rel- 
ativistic fluid. However, except near and within the cavi- 
ties, the total pressure P+P c (dotted lines) deviates very 
little from the pressure profile in the cluster before the 
cavity was introduced (dot-dashed lines). This is a con- 
sequence of the largely subsonic character of the flows. 
After about 0.1 Gyrs in both calculations the cosmic rays 
diffuse sufficiently so that the pressure ratio P c /P <C 0.2, 
the current threshold of 7-ray detectability of cosmic ray 
protons if nearby clusters are completely filled with cos- 
mic rays (Pfrommcr & Ensslin et al. 2004). Although 
the cosmic ray pressure and energy density are negligi- 
ble at these late times, we show below that cosmic rays 
continue to displace about the same volume of gas. This 
results in a very long-lasting global expansion of cluster 
gas that can be seen by the small discrepancies in the 
z-axis gas pressure profile (solid lines) relative to that in 
the initial cluster (dash-dotted lines) at times 0.3 Gyrs in 
Figure 3. Nevertheless, the close similarity of the initial 
and final gas pressure along the z-axis (and elsewhere) at 
late times is consistent with our finding that the gas den- 
sity and temperature (and therefore entropy 1 ) gradients 
are also only slightly affected by the cavity evolution. 

4. ENERGETICS AND EVOLUTION 
OF X-RAY CAVITIES 

Figure 4 shows the evolution of cluster gas energies re- 
sulting from cavity-lobe formation with cosmic rays. The 
four panels show the evolution as the cosmic ray diffu- 
sion coefficient k decreases with n e0 (in cm~ 3 ) over a 
wide range: n e o = 6x 10~ 3 (panel a), 6 x 10~ 4 (panel 
b), 6 x 10" 5 (panel c), and 6 x 10~ 6 (panel d). The total 
cosmic ray energy integrated over the entire computa- 
tional region E cr = J e c dV is shown with long dashed 
lines. E cr (t) is seen to rise until t cav = 2 x 10 7 yrs as 
the cavities form, then remain approximately constant 
after the cosmic ray source is turned off. E cr is not a 
strictly conserved energy. The small decrease in E cr (t) 
visible at times t > t cav can occur if the cosmic ray en- 
ergy density is reduced by a secular advective expan- 
sion with the cluster gas. This decrease in E cr (t) after 
t C av is stronger when k is smaller (panels a — > d) since 
the cosmic rays are more confined near the z-axis where 
most of the gas expansion occurs. The (small) total ki- 

1 Figure 8 shows (for the evolution in Figure 2) an approxi- 
mate z-axis gas entropy profile at 0.9 Gyrs relative to the initial 
atmosphere. When plotted, s(z) at 0.9 Gyrs is almost identical 
to that in the original atmosphere from the origin to z sa 25 kpc 
and beyond z ~ 42 kpc. In the intermediate region 25 < z < 42 
kpc the entropy at 0.9 Gyrs is depressed below the original profile 
approximately as shown in Figure 8. 
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netic energy in the cluster gas E^ in = 0.5 J p(u 2 +u 2 )dV 
is shown with dotted lines in each panel. The change 
in potential energy compared to the initial Virgo atmo- 
sphere, AEp t — J 4>[p{t = 0) — p]dV is shown with 
short dashed lines. The gravitational potential 4>(R) is 
found from the initial M87/ Virgo atmosphere by inte- 
grating d(j>/dR = -{l/p)dP/dR where R = (r 2 + z 2 ) 1 / 2 
is the radial coordinate. The change in the total gas ther- 
mal energy relative to the original atmosphere AE t h — 
J[e — e(t = 0)]dV is shown with (the lower) solid lines in 
each panel. Finally the dash-dotted line shows the total 
energy E tot = AE th + AE pot + E kin + E cr which is con- 
stant after time t — t cav to an excellent approximation 
and is equal to E cav = 10 58 ergs as expected. 

Also shown in each panel of Figure 4 is an approx- 
imate evaluation of the quantity 4PV where V is the 
volume of the X-ray cavity at any time (arbitrarily de- 
fined as the sum over all grid zone volumes in which 
the gas density is less than p(t = 0)/3) and P is an 
estimate of the average pressure in the (sometimes non- 
contiguous) zones containing the cavity. When cosmic 
rays are strongly trapped within the cavity, we expect 
4Py = E cav = E to t- In panel a of Figure 4 we see that 
4Py < E cav which can be expected for larger re when 
cosmic rays diffuse through the cavity walls. However, 
4Py > E cav is apparent in panels c and d, although APV 
never exceeds E cav = 10 58 ergs by more than about 20%. 
This may be due to our rather approximate estimate of 
P and V or it may be a real inertial overshoot just after 
the cavity is formed. In any case, we include 4PV(i) in 
Figure 4 because this is the cluster gas heating energy 
proposed by Rafferty et al. (2006) and McNamara & 
Nulsen (2007) in their discussions of cavity energetics. 

Figure 5 shows the approximate evolution of the cav- 
ity radius r cav (t) and its mean radius in the cluster 
Rcav(t). The cavity radius is found from the estimated 
cavity volume V by assuming that the cavity is spherical, 
r C av = (SV/An) 1 / 3 . The four lines for r cav and R cav in 
Figure 5 correspond to the four cosmic ray diffusivities 
re as it decreases in each panel in Figure 4: a, dotted 
line; 6, short dashed line; c, long dashed line; d, solid 
line. The buoyant trajectory of the cavities R ca v(t) is 
independent of re, but the cavities progress further along 
Rcav{t) when re is smaller. The cavity lifetime is deter- 
mined when r cav — ► 0. Notice that the cavity lifetimes 
are similar for n e o = 6 x 10 (long dashed line), and 
n e o = 6 x 10~ 6 (solid line), suggesting that the cavity 
evolution would not change much if n e o (and therefore 
re) were reduced further. 

The main result we wish to emphasize is the evolution 
of the cluster gas thermal energy AE t h(t). At early times 
during cavity formation when t < t cavi AE t h increases 
because of heating due to the weak shock that propa- 
gates away from the expanding cavity. This initial shock 
heating increases as re decreases (panels a — ► d in Figure 
4) because (for fixed E cav and t cav ) the shock strength 
increases when cosmic rays are more confined within the 
cavity. This is consistent with the ID cavity evolution 
described by Mathews & Brighenti (2007a). However, 
after the cavity is formed (t > t cav ) AE t h(t) decreases 
and becomes negative after t « 3 — 5 x 10 7 yrs. We 
see in Figure 4 that the final energy separation between 
AE t h and AE pot is independent of re. Also notice that 



in each panel of Figure 4 the average final value of the 
energy change 0.b(AE t h + AE pot ) is not zero, but is ap- 
proximately equal to the peak thermal energy acquired 
from the shock, ~ AE t h(t cav )- Although the total en- 
ergy change in the cluster AE pot + AE t h +Eki n increases 
as re decreases, the cluster gas experiences a net cooling 
(AEth < 0) when cavities are formed with cosmic rays. 
This global cooling is exactly opposite to the results dis- 
cussed by other authors in which cavity formation is re- 
garded as an important source of AGN heating required 
to balance radiative losses in the cluster gas and reduce 
central cooling in cooling flows. 

The energy evolutions shown in Figure 4 are not qual- 
itatively altered when the total cavity energy E cav is in- 
creased. For example when E cav = 10 59 ergs (deposited 
at (r, z) — (0, 10) kpc in time t cav = 2 x 10 7 yrs with 
n e o = 6 x 10~ 6 cm~ 3 ) all energies are larger but the 
proportions are similar to those in Figure 4. Although 
the shock is stronger in this case, so is the buoyant cool- 
ing and AEth still becomes negative after ~ 10 8 yrs. 
The maximum equivalent spherical radius of this high 
energy cavity is about r cav = 15 kpc at logt ~ 8.2, but 
the cavity disappears soon afterward at logt = 8.4 when 
Rcav = 55 kpc. More energetic cavities are bigger and 
more buoyant, but don't last proportionally longer. 

5. GLOBAL EXPANSION AND MASS OUTFLOW 

To verify that the cluster gas has in fact undergone a 
net expansion as a result of cavity formation, we com- 
puted the evolution of the cumulative gas mass distri- 
bution M(R). First we sorted the computational zones 
in cylindrical coordinates in order of increasing spheri- 
cal radius R = (r 2 + z 2 ) 1 ! 2 then we integrated over the 
7855 sorted zones within R — 50 kpc to determine M(R). 
Figure 6 shows M(R) in (one hemisphere of) the initial 
cluster (t = 0) and after t = 0.9 Gyr for two limiting 
values of re (n e0 = 6 x 10~ 3 and 6 x 10~ 6 cm -3 ). To 
more clearly illustrate the evolution of M(R) most of 
the R- variation in Figure 6 has been removed by plotting 
M(R)(R/10 kpc) -2 . The noise in these plots of M(R) 
arises because of the finite number of computational grid 
zones. Figure 6 clearly shows that mass has been re- 
moved from the inner regions of the cluster as a result of 
cavity formation, i.e. the cluster gas has experienced a 
net expansion. This expansion is remarkably insensitive 
to the cosmic ray diffusivity re, particularly at R 25 
kpc. The evolution of M(R) is expected to depend both 
on the cavity energy E cav and the location of the initial 
cavity in the cluster gas (10 kpc in our case). 

The largest radius plotted in Figure 6, R — 50 kpc, 
is essentially the cooling radius for the cluster, i.e. the 
radius where the local cooling time is 7.4 Gyrs, which 
is comparable to the cluster age. From this plot we 
estimate that AM W = 1.14 x 10 8 M© has been re- 
moved from within 10 kpc for both values of n e Q. At 
R = 50 kpc the mass outflow is AM 50 = 3.6 x 10 8 M© 
(for n e0 =6x 10~ 3 ) and AM 50 = 5.6 x 10 s M© (for 
n e0 = 6 x 10~ 6 ). Values of AM refer to the total cluster 
mass flow in both hemispheres, assuming that a pair of 
cavities are produced by symmetric double jets. 

Using the Ghizzardi et al. (2004) density and temper- 
ature profiles for Virgo, we estimate a feedback-free gas 
cooling rate of M c f « 8M©/yr. Therefore if cavities sim- 
ilar to the one we calculate here are formed with cosmic 
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rays every 6 x 10 7 yrs, this could result in a mass out- 
flow at R = 50 kpc equal to the cooling inflow at R = 
from radiation losses. When E cav = 10 59 ergs the total 
mass flowing beyond 50 kpc after time 0.9 Gyr is much 
larger, AM 50 = 4.6 x 10 9 M Q . (Note that the buoyant 
mass outflow scales with E cav , but the cavity lifetime 
does not.) Successive cavities in Virgo with E cav = 10 59 
erg appearing every 6 x 10 8 years could balance the mass 
flow from radiative cooling M c f w 8M Q /yr and also be 
consistent with the current absence of large cavities in 
M87/Virgo. Of course an expansion outflow at the cool- 
ing radius comparable to the steady state cooling rate 
M c f does not in itself shut down radiative cooling near 
the cluster center, but the cooling radius will become 
larger and the cooling time at every radius will increase 
due to the slightly lower gas density resulting from the 
global expansion. 

6. WHAT INCREASES THE POTENTIAL ENERGY 
AND DRIVES THE MASS OUTFLOW? 

The post-cavity mass flow has two aspects: (1) an in- 
crease in AEp t in Figure 4 and non-zero AM50 in Figure 
6 that result from the global lifting that occurs when a 
new cavity is formed, and (2) a mass circulation driven 
near the z- axis by buoyant gas containing cosmic rays. 
We discuss each in turn. 

6.1. Global Lifting 

Although gas certainly flows beyond R = 50 kpc in 
all our cavity evolutions, the fraction of total cosmic ray 
energy within this radius E cr ^/E cr varies. For example, 
at t = 0.9 Gyr and large k (n e o = 6 x 10~ 3 cm -3 ), 
the cosmic ray energy in R < 50 kpc, E cr ^o — 66.3 x 
10 56 ergs, is less than the total cosmic ray energy on the 
computational grid, E cr = 81.5 x 10 56 ergs. In this case 
about 20% of the cosmic rays have moved beyond 50 kpc. 
For smaller k (n e o = 6 x 10~ 3 cm -6 ) after this same time 
we find that -E cr ,5o = 41.8 x 10 56 ergs is identical to the 
total energy E cr , i.e. almost all cosmic rays are confined 
within 50 kpc. (Both E cr ^ and E cr refer to a single 
hemisphere.) However, the mass AM 50 flowing across 
R = 50 kpc is positive for both cases involving large and 
small diffusivities k. 

Evidently mass can flow across the cooling radius even 
in the absence of local buoyant transport. This must 
arise because of the quasi-adiabatic expansion of the en- 
tire cluster gas when cavities are formed. To understand 
this better, we begin by showing that the volume of gas 
displaced by cosmic rays V c depends on the local pres- 
sure P a in the cluster atmosphere and the total energy 
E cr of the cosmic rays, but not on the ratio of cosmic 
ray to gas partial pressures, P c /P, which changes as the 
cosmic rays diffuse into the gas. Imagine a uniform gas 
containing a spherical region of volume V that contains 
a total cosmic ray energy E cr . In pressure equilibrium 
P a = P c + P where P c = (j c - l)E cr /V and therefore 
V = (j c — l)E cr /P c . But the partial volume V c contain- 
ing gas that is actually displaced by the cosmic rays is 
V c = (P c /P a )V = (7 C - l)E cr /P a . Therefore as cosmic 
rays diffuse into larger volumes V, the total volume of gas 
that is displaced V c remains unchanged and independent 
of k. In a non-uniform cluster environment the volume 
of gas displaced by a given total energy of cosmic rays 
E cr varies with cluster radius as P a (R)~ 1 . 



AM50 can be estimated by assuming that all cosmic 
rays in the original cavity do not diffuse beyond their 
original cavity volume V cav and remain at their initial 
cluster radius R cav = 10 kpc. In this limit, and after 
several cluster sound crossing times, we expect that the 
mass flow across R = 50 kpc will be equal to the gas mass 
displaced by the original cavity, AM50 « V cav p{R cav ) = 
6.2 x 10 s M where V cav = (47r/3)r 3 a „, r cav « 6 kpc 
and p(R C av) = 4.72 x 10~ 26 gm s _1 is the original cluster 
density at R cav . This simple estimate of the outflowing 
mass is very similar to that found from Figure 6, AM 50 = 
5.6 x 10 8 M Q . 

However, discrepancies between computed and esti- 
mated masses AM50 are expected because of several ad- 
ditional details. The cosmic ray energy inside the esti- 
mated cavity volume is E cr — 3PioV cav — 101 x 10 56 
ergs, where P w = P a (R = 10 kpc) = 1.28 x lO" 10 dyne 
cm~ 2 . However, this energy is about twice as large as the 
value of E cr (t cav ) in the 1ow-k evolution shown in panel 
d of Figure 4, so the estimated value of AM 50 above 
should be reduced by a factor of ~ 2. However, AM 50 is 
expected to differ from V cav p(R ca v) because of changes 
in the mean volume V c of gas displaced by cosmic rays 
after their outward motion during time 0.9 Gyr. To es- 
timate this suppose that the e c -weighted mean radius of 
cosmic rays at time 0.9 Gyr is R = 30 kpc (see Figure 
2) where the gas pressure is P30 = 0.608 x 10~ 10 dyne 
cm~ 2 . For constant total cosmic ray energy we expect 
the volume of gas displaced to vary as V c oc P a ~ . There- 
fore if most of the cosmic rays are transported to ~ 30 
kpc, our original estimate of AM50 must be increased by 
~ P10/P30 = 2.1. The combination of both corrections 
leaves our original estimate nearly unchanged. Consider- 
ing the uncertainties in these estimates, we conclude that 
the total mass flow past R = 50 kpc found from Figure 
6 for the low-K evolution, AM 50 = 5.6 x 10 8 M , is 
consistent with a quasi-adiabatic expansion of the Virgo 
atmosphere past radius 50 kpc due to gas displaced by 
cosmic rays when the cavity formed and no buoyant mass 
flow across R = 50 kpc is required. 

6.2. Buoyant Gas Circulation 

In addition to this atmospheric lifting, buoyant mass 
flow largely within the cooling radius also contributes 
to the post-cavity cluster evolution. Buoyant gas in- 
creases AE po t and as it expands we expect AE t h to de- 
crease. But every buoyant element must be accompanied 
elsewhere in the hot gas atmosphere by a downward re- 
placement flow, largely balancing the buoyant changes in 
AE po t and AE t h ■ If cavities are created by feedback fol- 
lowing a small accretion onto the central AGN, a buoyant 
mass circulation flow is a natural response to the inflow 
expected from radiative cooling. The low-/t cavity evo- 
lution (n e0 = 6 x 10~ 6 cm~ 3 ) in Figure 2 is particularly 
interesting because all buoyant motion occurs within the 
cooling radius at 50 kpc. 

Figure 7 shows contours of the radial velocity v r (r,z) 
at four times during the low-/t cavity evolution in Figure 
2 (n e o = 6 x 10~ 6 cm -3 ). (A plot of the radial velocity in 
the high-K cavity evolution shown in Figure 1 is almost 
identical.) We plot v r (r, z) only near the z-axis since the 
radial velocities elsewhere in the atmosphere are very 
much smaller. Figure 7 shows a strong and sustained 
subsonic outflow until at least 0.1 Gyr. At t = 0.3 Gyr 
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a backflow occurs near the z-axis as the densest parts 
of the thermal "cavity jet" fall back toward the cluster 
center (see Mathews & Brighenti 2007b for more details). 
Later at t = 0.5 Gyr a mixture of positive and negative 
radial velocities is visible but all velocities are smaller 
and this trend toward quiescence continues until t = 1 
Gyr. 

Because the entropy decreases monotonically with ra- 
dius in the initial cluster atmosphere, we can use the 
entropy as a tracer to confirm that there has been a net 
mass outflow within 50 kpc during the 1ow-k evolution 
(which is essentially adiabatic apart from the weak shock 
wave). We define the gas entropy as s = 10~ 33 (e/p 5//3 ) 
in cgs units. In addition we expect that the buoyancy 
and outward flow of each gas element will increase with 
the partial pressure of the cosmic rays, P C /(P + P c ). In 
the upper panel of Figure 8 we plot P C /(P + P c ) against 
As = s(t = 0.9) — s atmi the change between the entropy 
s at 0.9 Gyr and the entropy at the same location in the 
original Virgo atmosphere, s atm . Each point corresponds 
to a computational grid zone. In those regions where 
P C /{P + Pc) is largest we see that As decreases system- 
atically with P C /(P + P c ), indicating that gas containing 
more cosmic rays at time 0.9 Gyr originally came from 
high entropy regions closer to the cluster center. 

The lower panel of Figure 8 shows (with small points) 
the location in the cluster at time 0.9 Gyr that contains 
gas with P C /(P + P c ) > 0.3 and the open squares mark 
those (significantly overlapping) zones where As < —0.1 
for the low-K cavity evolution. Gas with the highest cos- 
mic ray content is also the gas with the lowest entropy 
compared to that at the same radius in the original clus- 
ter. The spatial distribution of this low entropy gas is 
similar to the cosmic ray contours in the final panel of 
Figure 2. Low entropy regions in the lower panel of Fig- 
ure 8 also have a slightly smaller gas density than ad- 
jacent gas at the same cluster radius R without cosmic 
rays. During times of interest for the cluster evolution, 
the gas in these post-buoyant regions containing cosmic 
rays will never return to the cluster core where it initially 
resided. Figure 8 shows conclusively that low entropy gas 
originally at small cluster radii has been buoyantly trans- 
ported outward during the low-/t post-cavity evolution, 
but very little has moved beyond 50 kpc by time 0.9 Gyr. 
Cosmic rays from cavities are driving a mass circulation 
in the cluster gas out to a large fraction of the cooling 
radius. 

Evidently the long-term decrease in cluster gas mass 
AMio = 1.14 x 10 8 M Q within R = 10 kpc where the 
cavities formed arises due to the combined effects of the 
jet of thermal gas along the z-axis and the cosmic ray 
buoyancy that prevents much of this gas from return- 
ing back to R = 10 kpc. Recall that AM W is nearly 
independent of the cosmic ray diffusion coefficient. The 
cosmic ray energy required to remove a solar mass of gas 
from R = 10 kpc is about e© = E cav /AM W = 8.8 x 10 49 
erg Mq~ . If MSt = 8 M Q would cool each year in 
Virgo due to radiation losses alone, the cosmic ray lu- 
minosity required for a buoyant outflow to balance M is 
L cr = e Q M = 7.0 x 10 50 erg yr" 1 or 2.2 x 10 43 erg s" 1 . 
However, the rate that energy is created by accretion 
onto the - 3 x 10 9 M Q central black hole in M87/Virgo 
is L acc w O.lMc 2 = 4.6xl0 46 ergs~ 1 . The cosmic ray lu- 



minosity L cr required to balance a radiating cooling mass 
inflow in Virgo is ~ 2000 times less than the accretion 
luminosity generated by the black hole L acc . Therefore, 
if only a fraction e cr « 5 x 10~ 4 of the accretion energy 
converts to cosmic rays that inflate X-ray cavities, this 
might be sufficient to shut down the cooling flow. A sub- 
sequent calculation with multiple cavities and radiative 
cooling will be required to verify this. 

7. DISCUSSION 

In previous numerical simulations of X-ray cavity evo- 
lution (e.g. Bruggen & Kaiser 2002; Reynolds et al. 
2005; Gardini 2007; Pavlovski et al. 2007) as well as 
in the recent review of McNamara & Nulsen (2007) it 
has been concluded that the cavities have an important 
net heating effect on the cluster gas. This result can be 
understood because the simulated cavities were inflated 
with ultra-hot (but non-relativistic) gas. As these cav- 
ities are initially inflated, they perform an amount of 
heating ~ PV on the surrounding cluster gas, but the 
thermal energy inside the cavities also contributes to the 
cluster thermal energy budget. When we perform calcu- 
lations in which cavities are inflated with ultra-hot gas, 
we also find a (small) net heating when the total cluster 
thermal energy includes the energy injected inside the 
cavity. In this case the total energy introduced by the 
cavity is [7/(7 — 1)]-P^ = (2.5 — 4)PV, depending on the 
assumed value of 7 inside the cavity. (After times > 10 8 
yr most of the energy of hot-gas cavities is stored in the 
potential energy of the cluster.) The approximate values 
of APV plotted in Figure 4 show the estimated cavity 
heating when cavities are formed with ultra-hot gas (e.g. 
Birzan et al. 2004; Rafferty et al. 2007; McNamara & 
Nulsen 2007). 

But ultra-hot gas is not an appropriate substitute for 
cosmic rays. When the diffusion k is sufficiently small 
and cosmic rays are strongly trapped within the cavi- 
ties, we find here that the expanding cavities also heat 
the gas by ~ PV. This can be seen in Figure 4 where 
AE t h(t cav ) w PV. But the cosmic ray contents of our 
cavities never contribute to the cluster thermal energy. 
Instead, the displacement of cluster gas by the diffusing, 
buoyant cosmic rays results in an overall expansion and 
cooling of the cluster gas. 

It should also be recognized that the introduction of 
ultra-hot gas in previous numerical simulations of cavity 
evolution implicitly suggests that mass as well as cosmic 
ray energy is transported out from the central AGN, i.e. 
this conventional means of cavity formation must be re- 
garded as a (mass) circulation flow similar to those we 
have considered (Mathews et al. 2003; 2004). This as- 
sumed non-relativistic mass outflow can be quite large as 
we discussed in the Introduction. While we assume here 
that the AGN produces jets of pure relativistic energy, it 
is possible that real jets do in fact carry non-relativistic 
mass away from the central AGN (as in Brighenti & 
Mathews 2006). Nevertheless, current X-ray observa- 
tions have been unable to detect thermal emission from 
ultra-hot gas in the cavities (e.g. McNamara & Nulsen 
2007, but see Mazzotta et al. 2002). 

8. CONCLUSIONS 

It is generally agreed that the dominant source of feed- 
back energy in galaxy clusters is accretion onto central 
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black holes (AGNs). Some of this accretion energy is 
thought to be transported out into the cluster gas by 
jets, forming X-ray cavities. In previous numerical simu- 
lations of X-ray cavity evolution the cavities were formed 
by introducing ultra- hot, non-relativistic gas at some ra- 
dius in the cluster gas. 

For the first time we consider cavities that are formed 
exclusively by cosmic rays that can both diffuse through 
the cluster gas and be advected by it. We discuss the en- 
ergetics of X-ray cavity and radio lobe evolution in the 
hot gas within the Virgo galaxy cluster. For simplicity we 
consider cavities that are formed in time t cav = 2 x 10 7 
yrs by cosmic rays with total energy E cav = 10 58 ergs 
but which have a wide range of cosmic ray diffusion co- 
efficients k since this is currently the most uncertain pa- 
rameter in cavity-lobe evolution. For simplicity, we do 
not consider here secular energy losses in the cosmic rays 
due to synchrotron emission, Coulomb heating, inverse 
Compton losses, spallation etc., so the total cosmic ray 
energy can only change due to local expansion or com- 
pression by the thermal gas. In this approximation cos- 
mic ray pressure gradients interact with the cluster gas 
in a perfectly adiabatic fashion. 

We conclude that: 

(1) As cavities expand during formation, they generate 
weak shock waves that propagate into the surrounding 
cluster gas where the shock energy is dissipated. When 
the cosmic ray diffusion n is low, the maximum shock 
heating ~ PV ~ E cav /A is approached. As k increases 
the cosmic ray pressure gradients are lower and the shock 
heating is considerably less. In general, APV , which 
varies with time, is an inaccurate measure of E cav . These 
results are consistent with our earlier studies of ID cavity 
evolution in a uniform gas (Mathews & Brighcnti 2007a) . 

(2) The longevity of visible X-ray cavities and the radial 
distance that they move out in the cluster gas during 
their lifetimes both increase as the cosmic ray diffusion 
k decreases. Unlike cavities formed with hot gas (e.g. 
Gardini 2007), cosmic ray cavities remain coherent for 
> 10 8 yrs. 

(3) In spite of shock heating, cavities formed by cos- 
mic rays have a net cooling effect on the cluster gas. 
This is unlike the cluster heating when cavities are in- 
flated largely with ultra-hot but non-relativistic gas, as 
commonly assumed, that later contributes to the cluster 
thermal energy. 

(4) The qualitative character of post-cavity energetics 
remains similar when the cosmic ray diffusion coefficient 
k varies over a wide range. After several 10 8 yrs the 
total energy change in the cluster AE pot + AE t h + E^ in 
increases as n decreases. 

(5) As cosmic rays displace cluster gas, the entire clus- 



ter atmosphere expands and cools, increasing the cluster 
potential energy and decreasing its thermal energy. Suc- 
cessive cavities will result in an irreversible expansion of 
the cluster gas as long as most of the cosmic rays remain 
trapped within or near the cooling radius, as in our com- 
puted post-cavity evolutions. 

(6) Inhomogcneously distributed cosmic rays are an im- 
portant source of buoyancy and mass flow in the cluster 
gas. As the cluster gas surrounding cavities is penetrated 
by diffusing cosmic rays, this gas rises in the cluster po- 
tential, driving a net mass outflow in the cluster gas. 
Radial mass circulation occurs even when no gas flows 
out in the AGN jets. 

(7) After buoyant gas containing cosmic rays flows from 
the cluster core into distant, low-density regions, it re- 
mains there never to return during times comparable to 
the cluster age. Because of the necessarily positive en- 
tropy gradient in the pre-cavity cluster, the entropy in 
the post-buoyant gas in these distant regions is slightly 
lower than elsewhere at the same cluster radius. 

(8) To quantify the previous conclusion, suppose a mass 
AMj; within the cluster radius R where the X-ray cav- 
ity forms is irrevocably transported by cosmic ray buoy- 
ancy to distant regions in the cluster gas. We find that 
a cosmic ray energy of = E cav / AMr w 10 50 ergs 
removes a solar mass of gas from the cavity site. For 
example in the Virgo cluster a cosmic ray luminosity of 
only L cr = e Q M = 2 x 10 43 erg s _1 can remove gas at the 
cooling flow rate M w 8 Mq yr _1 that would otherwise 
occur in the absence of feedback energy. Only a fraction 
e cr w 5 x 10~ 4 of the accretion energy onto the central 
black hole in M87/Virgo, L acc w O.lMc 2 , if converted 
to cosmic rays that inflate X-ray cavities, is in principle 
sufficient to shut down the cooling flow. 

Our conclusions suggest that cosmic ray buoyancy pro- 
vides an important new means of understanding or pos- 
sibly solving the cooling flow problem with an outward 
mass circulation that may balance the inexorable inflow 
due to radiation losses. This outward mass flow due 
to cosmic ray buoyancy has the added attraction of not 
heating the gas too much - in fact cavity formation with 
cosmic rays produces a net cooling. 

It remains to be determined if the radio emission or 
7T° gamma emission from galaxy clusters are at levels 
that are consistent with the buoyant activity required to 
arrest the cooling flow rate. After a few 10 8 years most 
of the cosmic rays will be in very low density gas where 
such emission will be greatly reduced. 

Studies of the evolution of hot cluster gas at UC Santa 
Cruz are supported by NASA and NSF grants for which 
we are very grateful. 
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Fig. 1. — Evolution of an X-ray cavity in Virgo with cosmic ray diffusion parameter n e o = 6 X 10 3 cm 3 at six times shown in Gyrs 
at the upper right of each panel. Solid lines show the gas density contours p(r, z) in units of 10 — 26 g cm -3 . Dotted lines show with six 
contours the cosmic ray energy density e c (r, z) in units of 10 -12 erg cm -3 . Two adjacent contours are labeled and others can be found by 
extending the same additive variation. 
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Fig. 2. — Evolution of an X-ray cavity in Virgo with cosmic ray diffusion parameter n e o = 6 X 10 6 cm 3 at six times shown in Gyrs 
at the upper right of each panel. Solid lines show the gas density contours p(r, z) in units of 10 — 26 g cm -3 . Dotted lines show with six 
contours the cosmic ray energy density e c (r, z) in units of 10~ 12 erg cm -3 . Two adjacent contours are labeled and others can be found by 
extending the same additive variation. Due to crowding at early times the cosmic ray contours are difficult to distinguish in the first three 
panels. For the final three panels the outer two cosmic ray contours are 14.29 k, 31.43 (panel 0.300), 10.0 & 20.0 (panel 0.500), and 4.29 
& 11.43 (panel 0.900). 
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Fig. 3. — Pressure profiles along the z-axis in units of 10 -10 dynes cm -2 . Each panel contains four (often overlapping) profiles: the 
gas pressure P(r) (solid lines), cosmic ray pressure P c (r) (long dashed lines), the total pressure P + P c (dotted lines), and the initial gas 
pressure in the cluster before the cavity is introduced (dash-dotted lines). Upper four panels: show pressure profiles at four times for the 
evolution in Fig. 1 with n e Q = 6 X 10~ 3 cm~ 3 . Lower four panels: show pressure profiles at four times for the evolution in Fig. 2 with 
n e o = 6 X 10~ 6 cm~ 3 . Each panel is labeled with the time in Gyrs. 
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Fig. 4. — Global energy evolution in four cavities with decreasing cosmic ray diffusion coefficient k oc n e o characterized by density 
parameters n e o: 6 X 10~ 3 (panel a), 6 X 10~ 4 (panel b), 6 X 10~ 5 (panel c), and 6 X 10 -6 (panel d). The energies are labeled as follows: 
cosmic ray energy E cr (long dashed lines); change in potential energy AE po t (short dashed lines); kinetic energy E kin (dotted lines); 
change in thermal energy AE t h (lower solid lines); and the total energy Etot (dash-dotted lines). The total energy associated with the 
approximate cavity volume 4PV is shown in the upper solid lines. All energies are in units of 10 56 ergs and are those in the hemisphere 
containing our computational grid. 
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Fig. 5. — Approximate evolution of the cavity radius r cav (t) and its mean radius in the cluster R ca v(t). The cavity radius is found 
from the estimated cavity volume V by assuming that the cavity is spherical, r cav = [3V / 4,-k] 1 ' 3 . The four lines correspond to the four 
decreasing cosmic ray diffusivities K(n e ,n e o) in each panel of Figure 4: a, dotted line; b, short dashed line; c, long dashed line; d, solid line. 
Each overlapping trajectory R C av(t) ends at the time when r ca v — > 0. 
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Fig. 6. — Variation of the cumulative spherical mass gas distribution M(R) in one hemisphere of the cluster at times t = and 0.9 
Gyr for two values of the cosmic ray diffusion parameter n e o shown in parentheses. The cumulative mass is multiplied by (R/10 kpc) -2 to 
remove most of the radial variation between R = 6.3 and 50 kpc, the approximate Virgo cooling radius. 



16 




17 



Fig. 8. — Top panel: Correlation of the cosmic ray partial pressure P C /(P + P c ) with the entropy change As in each computational 
zone for the 1ow-k post-cavity flow over time t = 0.9 Gyrs. Each point represents a computational zone within cluster radius R = 50 kpc. 
Bottom panel: Location of regions in the cluster with high cosmic ray partial pressures (P C /(P + Pc) > 0.3) (points) and regions of large 
negative entropy change (As < —0.1) (open squares). 
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